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Abstract 

We evaluate the vector and scalar form factor of the pion in the presence of virtual 
photons at next-to-leading order in two-flavor chiral perturbation theory. We also 
consider the scalar and tensor pion form factors of the energy-momentum tensor. We 
find that the intrinsic electromagnetic corrections are of the expected size for the vector 
form factor and very small for the charged pion scalar form factor. Detector resolu- 
tion independent photon corrections reduce the vector radius by about one percent. 
The scalar radius of the neutral pion is reduced by two percent by isospin-breaking 
contributions. We perform infrared regularization by considering electron-positron an- 
nihilation into pions and the decay of a light Higgs boson into a pion pair. We discuss 
the detector resolution dependent contributions to the various form factors and pion 
radii. 
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1 Introduction 



In their seminal paper of 1984, Gasser and Leutwyler |I| set up the scheme to investigate the 
consequences of the broken chiral symmetry of QCD in a systematic fashion. This method 
is called chiral perturbation theory (CHPT) and rests on the fact that at low energies, the 
(pseudo)Goldstone bosons (the pions) interact weakly and saturate most n-point functions 
almost completely. The underlying power counting and an outline of the method had already 
been established by Weinberg in 1979 p[. In ref. [|l|, a whole series of applications was 
considered to one-loop accuracy (next-to-leading order), like e.g. the vector and scalar form 
factor of the pion, elastic pion-pion scattering, pion radiative decay and so on. By now, all of 
the observables and processes considered in ||l[] have been worked out to two-loop accuracy, 
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either by direct Feynman graph calculations or by making use of dispersive methods 
We mention here the pion mass and decay constant 0, the pion form factors ^, and 
TT7C scattering |Q [^l- In most cases these two-loop corrections are sizeably smaller than 
the corresponding one-loop terms, as expected by the power counting underlying CHPT. 
We remark that the improved two-loop representation of the vector form factor of the 
pion has been used to significantly decrease the hadronic uncertainty for the muon {g — 2) 
prediction In addition, there is a series of upcoming experiments to pin down precisely 
the TTTT S-wave scattering lengths (via Ki4^ decays at BNL and DA$NE or pionium lifetime 
measurements at CERN). However, as first pointed out in ref. and further quantified 



in ref. [0, electromagnetic corrections to elastic vrvr-scattering in the threshold region are 
of the same size as the hadronic two-loop contributions. At the level of the presently 
achieved accuracy, it is therefore mandatory to include virtual photons and work out the 
dual effects of isospin violation due to the light quark mass difference and the electromagnetic 
interaction. This machinery was essentially started in refs. |TT|, |T^, for the three flavor 



case in connection with the violation of Dashen's theorem (see also ref. [T^), which is of 
importance for extracting the quark mass ratios from the masses of the Goldstone bosons. 

Here, we will be concerned with the vector and scalar form factor of the pion as well as 
pionic matrix elements of the energy-momentum tensor in the presence of virtual photons. 
The scalar form factor is of interest for various reasons. Due to its quantum numbers, its 
phase is related to the strong tttt interaction in the isospin zero S-wave and it also exhibits 
the unitary cusp at the opening of the two-pion threshold. Furthermore, it gives a direct 
measure of the QCD quark mass term in the pion. In addition, the scalar form factor 
plays an essential role in the dispersive analysis of the so-called pion-nucleon cr-term, which 
allows one to extract information about the strangeness component in the nucleon wave 
function. In the pion-nucleon sector, isoscalar quantities can be very sensitive to strong 
and/or electromagnetic isospin breaking, see e.g. [1^, |T^, 0. It is therefore mandatory 



to evaluate the effects of virtual photons on this scalar-isoscalar quantity. Similarly, the 
hadronic two-loop representation of the vector form factor is extremely precise and has even 
been used to extract the pion charge radius from space- and time-like data at low energies . 

'^'^Note that using dispersive methods does in general not allow to study the quark mass dependence of 
observables. 
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Here, the celebrated AdemoUo-Gatto theorem [|T^ lets one expect that symmetry breaking 
corrections are fairly small, but their precise magnitude has so far not been worked out. 
Another set of form factors is related to the pion matrix elements of the energy-momentum 
tensor pO[, which play a role in the decay of a light Higgs into two pions 12^, ^ or heavy 



quarkonia transitions into the ground state with the emission of a pion pair. Of particular 
interest is the trace of the energy-momentum tensor which allows one to study the interplay 



of conformal and chiral symmetry |j23[ or the hadron mass composition [^]. Clearly, it is 
of interest to learn about the modifications induced by virtual photons, although we expect 
these to be small effects. 

Of course, whenever one deals with virtual photons, one encounters infrared (IR) singularities 
due to the vanishing photon mass. These divergences can be cured by taking into account the 
effects of radiation of soft photons which have an energy below a given detector resolution. 
For that, one has to consider cross sections or decay widths. In our case, we consider 
the squares of the corresponding form factors since that is the way they appear in the 
pertinent cross sections or decay rates. For the vector form factor, we consider the process 
e^e~ —>■ 7r+7r~. The less easily accessible scalar form factor together with the form factor of 
the trace anomaly can be studied in the decay of the Higgs boson into two pions. We note 
that to the order we are working, we have to consider only the radiation of exactly one soft 
photon from any one of the final-state pion legs. Radiation of two or more soft photons is 
of higher order in the fine structure constant. 

The paper is organized as follows. In sect. § we briefiy discuss the strong and electromagnetic 
effective chiral Lagrangians underlying our calculation. We work in SU (2) and to fourth order 
in the chiral expansion, counting the electric charge e as a small parameter like external 
momenta or pion mass insertions. In sect. ^ we give the results for the scalar and vector 
form factor of the pion at next-to-leading order including virtual photons and discuss the 
various checks on our calculation. The infrared regularization is performed by means of 
a finite photon mass. In sect. ^ we briefiy review the effective Lagrangian for pions and 
external fields coupled to gravity and include virtual photons. The pion matrix elements of 
the energy-momentum tensor are calculated in sect. § to fourth order in the chiral expansion 
including virtual photons and compared with previous results (obtained in the isospin limit). 
The problem of infrared regularization is taken up in sect. ^ Our results are presented and 
discussed in sect. 0. A short summary is given in sect. M. 



2 Effective Lagrangian with pions and virtual photons 

In this section we discuss the strong and electromagnetic effective Lagrangians underlying 
our calculations. The various pieces have already been discussed in the literature, so we will 
be brief and just collect the terms actually needed in what follows. The effective Lagrangian 
has the energy expansion 

Ces = C^^''''^ + C^^''''^ + + + . . . , (2.1) 
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where the superscripts 2, 4 refer to the chiral dimension and "str" and "em" denote the strong 
and the electromagnetic terms, respectively. The ellipsis stands for higher order terms not 
considered here. Our power counting is based on the observation that 6^/47? ~ M'^/ {AtiFY — 
1/100, with the pion mass and F the weak pion decay constant. This means that we 
consider the electric charge as a small parameter similar to the external momenta and pion 
mass insertions. Consequently, any term of the form e^" M^' with ra+m+Z = 2 is counted 
as fourth order. In what follows, we denote such a term as 0{q^) without differentiating 
between external momenta, pion mass insertions, or powers of the electric charge. 

At leading order, the strong Lagrangian is nothing but the non-linear a-model coupled to 
external sources, 

C^^^^^^) = !l^uX + X-,) , (2.2) 

where (. . .) denotes the trace in flavor space, U{x) = u'^{x) parametrizes the pion fields, 
Ufj^ = iu^D^Uu^ is the conventional pion viel-bein and the chiral covariant derivative 
in the presence of external vector and axial-vector fields. Furthermore, x = 2i?o(s + ip) 
subsumes the external scalar and pseudoscalar sources, and x+ = u^X^^ + ux^'^- The 
external scalar source contains the quark mass matrix, s{x) = diag(m„,m£;) + and 
Bq = 1(01^510)1/1^2 measures the strength of the spontaneous symmetry breaking. We as- 
sume -Bo ^ (standard CHPT). Throughout, we work in the so-called cr-model gauge, 
which simplifies the calculations. It is given by 



U{x) =(r{x)l2 + iT-n{x)/F , a{x) = ^l-7r^{x)/F^ , (2.3) 
with I2 the unit matrix in flavor space. The relevant terms of the strong Lagrangian at 



next-to-leading order read (we use the notation of ref. p5| ) 



= ^{X^r + '^{u,x'-) + ^(/rK,«.]) - ^(X-)^ , (2.4) 

with X- = u'^xu'' ~ ux'^u, X- = m^-D^xm"!" — uD^x^u, = uF^uU^ + F^^u, and F^^ = 
d^Ay — dvA^ is the conventional photon field strength tensor. The low energy constants 
(LECs) /s and /4 appear in the chiral expansion of the pion mass and pion decay constant 
(and thus in the scalar form factor), whereas features prominently in the pion charge 
radius as discussed below. Of particular interest to us is the term ~ /y since it embodies the 
strong isospin violation ~ {ma — rriu). After standard renormalization, the /j turn into scale- 
dependent, renormalized LECs called ^[(A). As first proposed in ref. |Q, we work instead 
with scale-independent LECs k obtained from the /[(A) via 



327r2 



k + loi 



A2 



(2.5) 



with A the scale of dimensional regularization and (3i the pertinent /^-functions. We remark 
that we use the neutral pion mass as the reference scale since the charged to neutral pion mass 
difference is almost entirely of electromagnetic origin. Note that 1^ does not get renormalized 
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{/3j = 0) and in that case we simply set li = IGtt^ Ij. This procedure is adopted for all other 
counterterms with vanishing /3-functions. The numerical values of the LECs will be given 
later. 



The leading order electromagnetic Lagrangian reads [T^ 

C^''''^^ = -\f^uF'"' -Ya^d-Af + C (QUQU^) , (2.6) 

with Q = diag(2/3, — 1/3) the quark charge matrix and a the gauge fixing parameter. In 
what follows, we work in the Feynman gauge a = 1.0 In the cr-model gauge for the pions, 
the last term of this Lagrangian leads only to a term quadratic in pion fields. Consequently, 
the LEG C can be calculated from the charged to neutral pion mass difference, 

(M'^^ - Ml^U = mlU = ^ = '^Ze'F" , Z=^^. (2.7) 

We have Z = 0.89 if we use for F its value in the chiral limit, Fq = 88MeV. At next-to- 
leading order, the relevant terms from the electromagnetic Lagrangian read 0, ^ 

£(4, em) _ kiF^{QUQU^y + k2F^QUQU^){D^UD^'U^) 

+ hF^ {{U^D^UQ)^ + {D^UU^Qf) + k^F^ {U^ D^UQ){D^UU^Q) 

+ krF^ {QUQU^){xU^ + Ux^) + k^F^ {{U\ - x^U){U^QUQ - QU^QU)) 

+ k^ F' (Q2) (QUQU^) + k^o F^ (Q2) i^D^UD^U^) + k^^ F^ (Q^) (;^f/t + f/^t ) 

+ k^^F^ixU^ + Ux^)Q + {x^U + U\)Q){Q) . (2.8) 

Note that the covariant derivative contains external vector sources and the virtual photon 
fields, see e.g. refs. |l^, 0. The terms proportional to fci and kg are only needed for the mass 
renormalization of the charged pions and the trace of the energy-momentum tensor (as 
discussed below). For estimates of some of the LECs ki needed in connection with the 
violation of Dashen's theorem, see refs. As shown in ref. []TU[, one needs additional 

divergent terms in the one-loop generating functional to renormalize the photon propagator, 

(2.9) 

(where e = d — 4) in terms of left /right-handed gauge fields (for more precise definitions, we 
refer to app. A of ref. |10[)- In our case, we can substitute G^i^ by QF^^ and set U{x) = 1. 
Thus, eq. (2^) boils down to one term of the form 

£(4,em,7) = F^^F^^- ^ (2.10) 

with the corresponding /3-function being = 2/3. From that, one easily deduces the 
Z-factor of the photon. 



^7 = 1+ 3(^7-1), (2.11) 

in terms of the renormalized, scale-independent low energy constant k.y. 

#'*Here, we have already set Ql = Qb. = Q since we are not concerned with the construction of the effective 
Lagrangian. 
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3 Pion form factors 



The pion vertex functions with one external vector or scalar source are parametrized in 
terms of the vector and scalar form factor of the pion, respectively. The vector form factor 
is experimentally accessible in electron scattering off pions or electron-positron annihilation 
into a pair of pions. Due to the quantum numbers, the first non-Goldstone meson, the 
p(770), features prominently in this channel and limits the range of applicability of CHPT. 
However, even at much lower energies, the tail of the p is visible, it saturates the LEG Iq 
and thus largely determines the pion charge radius, as discussed below. The scalar form 
factor can only be inferred indirectly due to the absence of any external scalar isoscalar 
source. It plays a central role in the discussion of low energy theorems since the scalar 
radius is intimately related to the change of the pion decay constant when the quark masses 
are turned on, as pointed out in ref. We now turn to a more detailed discussion of these 
two form factors. 



3.1 Vector form factor 

The vector form factor of the pion is defined by 

{7:\p')Ap)oni\V^^\Q) =ie^"^{p'^-p,)F:{s) , (3.1) 

with s = [p' + pY and the c-component of the quark isovector vector current, = 
g7^(r'^/2)g. Gauge invariance requires F^{0) = 1. Since the neutral pion is its own an- 
tiparticle, the corresponding form factor vanishes identically. At next-to-leading order in the 
chiral expansion, there is no strong isospin breaking in |jl9|. The low energy expansion 



of the vector form factor defines the electromagnetic charge radius of the pion, 

F^{s) = l + ^-{r')]:s + 0{s'). (3.2) 

The one-loop representation of F^ has been given in whereas the expansion to two loops 
was worked out using dispersive methods in 0] and by direct diagrammatic calculation in p . 
Here, we are concerned with the virtual photon effects at next-to-leading order (one-loop), 
including in particular virtual photon loops. 

The pertinent Feynman diagrams for a complete calculation at fourth order with virtual 
photons are depicted in fig. 0. In the first row, we show the strong one-loop graph together 
with the irreducible one-photon loop diagrams. The second row contains the reducible 
photon loop diagrams, below the standard wave function renormalization as well as strong 
and electromagnetic counterterms are depicted. The graphical representation of the Z-factor 
is also shown. We give here the Z-factors for the neutral and charged pions (note that Z^o 
is needed for the calculation of the scalar form factor in the next paragraph), 

Z^o = l-^ + e''(2{2h + h)-'^{k2 + ho)y (3.3) 
+ 2e'77^-2 — l + 21og— ^ --e2(A;2 + M , (3-4) 



F2 M2± V47ry V M^±J 9 
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z = 



Figure 1: Vector form factor of the pion to one loop. Solid and wiggly lines denote pions and 
photons, in order. Strong and electromagnetic counterterms are depicted by the filled box 
and circle, respectively. Diagrams that vanish in dimensional regularization are not shown. 

with 

A,.=2Mj.(L + ^log^^). (3.5) 

Here, L contains the poles in l/{d — 4) as (i — 4. Note that the fcj in eqs. ( p.3| , p.4|) are the 
non-renormalized LECs of the dimension four counterterms. The photon loop graphs lead to 
a complication, namely they are not infrared finite. To handle these IR divergences related 
to the zero photon mass, we introduce a small photon mass m^, giving a modified photon 
propagator 

Observables are, of course, independent of m^, as discussed in section ^. At fourth order, 
the vector form factor of the pion takes the form: 

= 1 + (4^ + 



with 



V47ry 13 ml±-s ' 4M2±-s 3V s ^^^ 6 

+ + ^ ^-^] l°gXr^ ' 3-7 

Mt± \ sa a — 1 J M^± J 



/ 4M2+ 
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and 



Ka{s) = £ dx\og(l ~ x{l - x)^y a = {0,±}, (3.9) 



1 log(l — x(l — x) 



TT 



G{s) = I dx^ • (3.10) 



1 — x(l — xl 



77=1= 



We note that the irreducible photon loop graphs are singular at s = 4M^± whereas the 
reducible photon loops are smoothly varying with energy. A representation of G{s) in terms 
of dilog-functions is given in ref. (since G{s) is related to the IR-finite part of their loop 
function G_|__^(s)). The non-renormalization theorem of the electric charge is fulfilled since 
all the various contributions are of 0{s). This is easily verified using 

In particular, the electromagnetic counterterms from the fourth order Lagrangian, which 
are momentum-independent, cancel exactly with the corresponding counterterms in the 
Z-factor, compare fig. ^ Concerning the strong contribution, it is well-known that the 
counterterm ~ Iq dominates the pion charge radius. This LEG is entirely saturated by the 
p-meson, which means that even in the theory of pions only, this resonance plays a visible 
role. 



3.2 Scalar form factor 

The scalar form factor of the pion is defined by 

{^''{p)^\p)out\muUU + mddd\^) = 5''''T^{s) . (3.12) 

In order to display the momentum dependence, one introduces the normalized scalar form 
factor r^(s) = Tt,{s)/Tt,{0) so that r^(0) = 1. The low energy expansion of the normalized 
scalar form factor defines the scalar radius of the pion, 

T^{s) = l+^-{rY^s + 0{s'). (3.13) 

Note that while the electromagnetic radius is directly measurable in electron scattering, the 
absence of scalar-isoscalar external sources allows only for an indirect determination of the 
scalar radius by e.g. a dispersive analysis of elastic tttt scattering |21]. At s = 0, the scalar 
form factor is proportional to the pion expectation value of the QCD quark mass term 

"^QCD = muUU + rriddd = -{rriu + md){uu + dd) + -{rriu — md){uu — dd) 

— m{uu + dd) + rh{uu — dd), (3-14) 
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in terms of the isoscalar and isovector components. Therefore one can apply the Feynman- 
Hellman theorem ^ separately to the isoscalar term, 



and similarly to the isovector one, 



ml 

dm 



dMl 
dm 



(7r|gg|7r) 



(7r|gr^g|7r) 



(3.15) 



(3.16) 



In the case of virtual photons, the arguments leading to eqs. ( |3.15| , |3.16| ) are not altered. 
These equations thus serve as a powerful check on our calculations. The one-loop represen- 
tation of has been given in ||1[, whereas the expansion to two loops was worked out using 
dispersive methods in [^] and by direct diagrammatic calculation in 0. We are concerned 
with the virtual photon effects at next-to-leading order. To the order we are working, strong 
isospin violation only occurs in r^(0). 



o 




Figure 2: Scalar form factor of the pion to one loop. The double line denotes the coupling 
to the scalar-isoscalar source. For further notations, see fig. |l[ 



The pertinent tree level and one-loop Feynman diagrams needed to calculate the scalar form 
factor at 0{q^) are shown in fig. ^. The scalar form factor for the charged pions is given by 



F^ifs) 



f.±(0) 



1 



1 



s I 



log 



(47rF)2 



1 



{Ml,-s)K^{s) 



+4 



2M\ , a+l 



a - 1 

2 

2 



S O 



(47rF)2 2 



-2Ze^F'^{ 1 + 



log 

-) 

47r/ 



1 ) log ""^ 



M'o (3 + 4Z) 



-V 

47ry 



4-(6 + 8Z) log-^ + ^ 



(3.17) 



(3.18) 
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and similarly for the scalar form factor of the neutral pions 



r.o(.) = l + ^^|s(r4-l-21og^) + (M,^o-.)ir±(.)-^Ko(.)} , 

(3.19) 



with 



k'^ = lfl + 2Z + 8Z')h-'^Z%-l(5 + Z + 4Z')h + ^(27 + 4Z)zho{?-2i: 



6V / " 9 6 

rriu + rrid 



-0.3 , (3.22) 



and the loop functions Ka{s), G{s) as given in the preceding paragraph. The light quark mass 
difference only appears in the normalization of the scalar form factor. This is related to the 
fact that isospin breaking only appears at next-to-leading order, see the last term in eq. (p.4| ). 
This operator only has non-vanishing matrix-elements with two neutral pions (similar to 
the operators C5 and /2 which appear in neutral pion scattering off nucleons |jl6|, |l^]). To 
arrive at the formulae for r^0 7r±, we have made use of the electromagnetic renormalization 
of the pion masses as given in refs. 0, [10 . 



4 Effective Lagrangian coupled to gravity 



The energy-momentum tensor for chiral effective theories has been studied in ref. [^. This 
extended Lagrangian can be constructed by including the metric tensor as an external source 
which couples to the energy-momentum tensor. At next-to-leading order, one has three new 
terms (plus additional terms without pions, needed for the renormalization). The energy- 
momentum tensor plays a role in the pionic decay channels of a hypothetical light Higgs 
particle via the matrix element (ttttIB^IO). We will later use this process to perform the IR 
regularization of the scalar pion form factor. Here, we briefly collect the pertinent results 
from ref. [^, adapted to our notation. In addition, we also include virtual photons. 

We work in SU (2) and rewrite the terms given in ref. pOI in terms of LECs /ii,i2,i3 (note 



that we keep the numbering as in [|20[). The pertinent strong Lagrangian coupled to gravity 



reads (terms just needed for renormalization are not shown )p^ 

^(4,str,R) ^ ^ {D^UD'^U^) + 1,2 R"" (D^UDM^) + h, R {xU^ + Ux^) , (4.1) 
where Rfj^i, and R are the Ricci tensor and curvature scalar, respectively, 

Rf,u = RnXu , R = Rn ^ (4.2) 



At leading order one has no new terms. One simply raises and lowers the Lorentz indices of the non-linear 
cr-model with the metric tensor y^j/. 
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derived from the curvature tensor and the Christoffel symbols, 



pA o -pA a -pA _|_ -pA -pa pA pa /a n\ 

= ^ g'^ {d,g,. + d,g,„ - d^g,,) . (4.4) 

From these formulae it becomes obvious why R and R^u count as 0{q^). For the case of 
SU{2), the /^-functions for In and /13 are 

Ai = ^, Pi3 = l- (4.5) 
8 

Note that the coupling /12 is not renormalized. For an estimate of the LECs and a more 
detailed discussion, see ref. ||2^. As before, we introduce scale-independent LECs rii,i2,i3 as 
described in sect. |^. As numerical values, we use 

/n = 8.7, ri2 = -0.4, [13 = 9.4. (4.6) 

These values have been obtained from the ones given in ref. by running them down to A = 
M^o using the SU{3) /5-functions and converting the obtained values for the renormalized, 
scale-dependent LECs into the scale-independent ones according to eq. (^.51). 



We are interested in the virtual photon corrections to the pionic matrix elements of the 
energy-momentum tensor. There is exactly one novel electromagnetic counterterm at fourth 
order, 

£{4,em,R) _ p2 p, ^QUQU^) . (4.7) 

Performing standard renormalization, the /3-function related to the LEC fcis is 

Pi, = ^ + lz. (4.8) 

It is obvious from the power counting that there is only one such term. The electric charge 
has to come in squared, and since R is of order g^, we cannot have derivatives on the pion 
fields or quark mass insertions. 



5 Pion form factors of the energy— momentum tensor 



Consider now matrix elements of the energy-momentum tensor in one-particle states like 
e.g. the pion. The general form of such matrix elements follows from Lorentz invariance and 



four-momentum conservation [20 



{n%p')\Q,,\n\p)) = 6'' 



^{sg^^ - q^q^) Qi{s) + ^P^^P^. 62(5) 



(5.1) 



with Pfj, = p'^+ p^, q^ = p'^ — Pfj,, and s = q^. The two form factors reflect the fact that the 
energy-momentum tensor contains a scalar and a tensor part. In what follows, we are mostly 
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interested in the matrix element of the trace of the energy-momentum tensor, G^, which is 
a scalar. The corresponding matrix element takes the form (neglecting isospin indices) 



(7r|0;:|7r) = 0.(.) = e,{s) + ^{4Ml - s) Q^is) . 
Because 0oo is the Hamiltonian density, we require the normalization 02(0) = 1. 



(5.2) 



O 



9 




9 9 
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Figure 3: Pion form factors of the energy-momentum tensor to one loop. The double line 
with the circle-cross denotes the coupling to the energy-momentum tensor. For further 
notations, see fig. 0. 

The pertinent Feynman diagrams needed to evaluate the form factors of the energy-momen- 
tum tensor at next-to-leading order (including virtual photons) in one-pion states are shown 
in fig. ^. At tree level, one simply has 0i(s) = 02(s) = 1. The full one-loop result for the 
charged pions reads 



1 - ™ {- - 2. log ^ + -(. + Ml.) Ms) + -(. + 2MI.) K^is) 



-ImIo ( ^Ko(s) + 1)- 2s(% + h. 



+ ( f V + Z)KAs) + 1{2Z - 1) ( + I 



S-2MI 



^ G{s) + 2Ml^ (h{s) - 2H,{s)) (5.3) 
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3 V M\ / V sa " a-l J " M- 
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efis) 



2s 
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4nJ 



4s 



+ 6- 



3^ 



/Cl5 + 2 log 



(5.4) 



4M2, 
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2 



2M2, 



G{s) 



log 



3s 
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log- 
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1 log 



(5.5) 



and similarly for the neutral pions we find 
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,(5.6) 
(5.7) 
(5.8) 



in terms of the loop functions Kq^±{s), G{s) defined in eqs. (|3.9| , p.lOj) and 



i/(s) = — 
s cr 

1 



HAs) 



a" 



with 



- log 

2 ^ 



7-± 



M2, 



lot 



(J 



\ogt_ 

l-t 



log 



M2, 



?7r 



l±a-2M2±/s 



1±(T 



(5.9) 
(5.10) 

(5.11) 



The following remarks are in order. First, in the isospin limit (m„ = m^, = 0) we recover 
the results of ref. pO[. The normalization of 62(0) is not affected by the photon loops, in 
particular not by the IR divergent terms. We also note that for the charged pions, Bi diverges 
logarithmically at s = 0. We have therefore only given the finite part of (0), denoted by 
the superscript "fin" . This divergence is due to the photon tadpole and self-energy graphs 
with coupling of the energy-momentum tensor to the photon line. Such diagrams are not 
present in the case of the scalar pion form factor simply because a scalar-isoscalar source 
has no coupling to virtual photons only. We remark that after mass renormalization, no 
counterterms from and contribute to Finally, we note that the new loop 

functions H{s) and Hi{s) cancel in the matrix elements of the trace of the energy-momentum 
tensor as defined in eq. ( |5.2|) . Consequently, these matrix elements are well behaved at s = 0. 
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We can now combine the results collected in eqs. (^^^) to work out the pionic form factors 



of the trace of the energy-momentum tensor. For the charged pions this leads to 

+ {\{Ml. - 2s)s + .(3. - 2Ml,) log ^ - - Ml,){s + 2Ml,) Ko{s) 

As{s + 2M2±) K±{s) + s{s - 2Ml,) In + As{s - M^o) 1^ + ^sM^^o } 
+ (|-)'{(2M,2,+s)(l + 4Z)i^±(s) 



+ (2M,2, + s f ^fj' Gis) + (3 + AZ)s (k,, + 2 log ^ 
+4 (2M.^. + .) C-^^ log ^ - 1 Vog , (5.12) 



s 



sa a -I J M,, 

and similarly for the neutral pions 

e^o(s) = 2Mlo + s 

+ - 2s)s - 2s{s - Ml,) log ^ - \mI,{s + 2Ml,) K,[ 

- {s- Mlo){s + 2M^„) + s{s - 2Mlo) hi + 4s(s - M^o) /"is + ^sM^o I13 

(5.13) 

One process in which these form factors play a role is the decay of a light Higgs boson into 
a pion pair. The corresponding matrix element takes the form 



{7T\p)7T\p')\C,s\H) = --g{s)6''' , g{s) = lQ^is) + l(t^is) + S^is)) , (5.14) 

where H denotes the Higgs field and the form factor S*^ is related to the expectation value of 
the strangeness operator rrisss in the pion. In what follows, we will neglect this form factor. 
Also, V is the scale of electroweak symmetry breaking, v~'^ = ^j2Gp ~ (250GeV)~^ (with 
Gp the Fermi constant). The width of the Higgs decaying into a pair of charged pions is 
then given by 



with niH the Higgs mass. For more details, we refer to ref. [21 . 



6 Infrared regularization 

So far we have dealt with the IR divergences by introducing a small photon mass m^. We are 
now going to cure this by considering processes with additional soft photon radiation from 
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the final-state pions with photon energies below a given detector resolution AE. Before 
discussing the details for the various form factors, we have to address the power counting. In 
fact, for the treatment of IR divergences, only the counting in e (or, equivalently, in the fine 
structure constant a) is of relevance. Since we consider one-photon loop diagrams, these 
scale as e^. Any radiation of a single soft photon carries a factor of e in a given amplitude, 
so that, for a pertinent cross section (which is the square of such an amplitude), we get 
contributions of order e^. Radiation of multiple soft photons is therefore related to higher 
order photon loop diagrams. The procedure we have to follow is standard, for a textbook 
treatment we refer e.g. to refs. []30| , |3T[| . Consequently, we do not concern ourselves here with 
initial state radiation. 



6.1 Vector form factor 




Figure 4: IR regularization of the vector form factor to one loop. The photon radiated from 
the charged pion lines is soft. 

Consider now the vector form factor. As shown in fig. ^, it can be measured e.g. in electron- 
positron annihilation. The corresponding tree graphs with soft photon radiation are also 
depicted in that figure. We evaluate the square of the vector form factor. After the IR 
regularization procedure, it takes the form 



e Y fl28M2± - 13s , , 4s 4: / MK _ _ . 1 



+2 — — —K±(s) -{-^K±(s) + -] (Q.i] 

\AnJ l3 4M2±-s 4M2±-s 3V s ^ ^ 6/ ^ ' 
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with the resolution function Ry{s, AE) 



Rv{s, AE) 



s — rat. 
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s + 2m1 s 



AE 



dll 



\ 
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4M2, 



- 12 



MK 1 



s + 2m1 s cr^ 



dlla{l) 



AE 



4mg s 



2M\ 



2m^ 



(i/ / log 



l + a{l) 
l-a(/) 



dl log 



_2/_ 

l + a(/) 
1 - a{l) 



4M\ /-AE dl 



s a 



a{l) 



s a 



+ 



2MK 



s a 



I Vl-a2(/) 

ci/ (l + a{l) 1 



log 



1 - a{l) l + a 



where 



a{l) 



s - 4M2± - 2^sl 



(6.2) 



(6.3) 



\ s - 2^sl 

cr(0) = cr as defined in eq. (|3.8|), and rrie is the electron mass. AE is the detector resolution 
which should be supplied by the experimenters. Typically, AE = 10 ... 20 MeV. The result- 
ing expression is obviously IR finite. For later discussion, we call the term ~ \og{AE / M^±) 
the IR finite contribution and the term proportional to the resolution function the finite 
bremsstrahlung contribution. Finally, we remark that if one chooses a different process to 
perform the IR regularization, the resolution function would be different. 



6.2 Higgs decay and scalar form factors 

We now turn to the IR regularization of the scalar form factors. Let us genuinely denote 
by St,±{s) any one of the charged pion scalar form factors r7r±(s), 0^±(s), or any linear 
combination thereof. To the order we are working, the regularization can be performed by 
considering a physical process like the decay of a light Higgs H into two pions accompanied 
by soft photon radiation from the final-state pions. The isospin-conserving part of this 



process has been studied in great detail in ref. see also ref. [|22|. We have used that 
calculation as a check whenever possible. Of course, the calculation is somewhat academic 
since such a light Higgs seems to be ruled out by experiment. It is, however, the simplest 
process which allows us to perform the IR regularization. Alternatively, one could consider 
processes like ip' —>■ tpTTir. More precisely, we only need to consider the emission of one soft 
photon from any one of the charged pions in the final-state. The method is depicted in fig. ^, 
which shows the tree graph for the H —>■ tttt decay via the corresponding scalar form factor 
and the two graphs with one soft photon. Denote by S^±{s) the leading term in the chiral 
expansion of any scalar form factor and by S^±{s) the next-to-leading order contribution. 
To one-loop accuracy, the IR-divergent part of the amplitude takes the form 

sg"'W = 4(^)'e'W«'og-;^. (6.4) 
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TC"" 71" 71"" 71" 71"" 71 




H H H 

Figure 5: IR regular ization of any scalar form factor ^^±(5) to one loop. The solid line 
labeled "H" denotes the light Higgs particle. The photon radiated from the charged pion 
lines is soft. 



with 

i= ^ log 7-1- (6.5) 

sa a — 1 

With that, we can write the square of the form factor as 

+2Re5i°2(s) S^X"\s) + 0{e^) , (6.6) 
with the resolution function Rs{s,AE) given by 

^ / A 16 r _2 f^^ dl / a(l) sa \ ^ /"^^ „ , /l + a(/)~ 

2 /-^^rf/ ,„„^l+^(0 



where a{l) has been defined in eq. ( p.3| ) and (j(0) = a in eq. (p.8|), respectively. Again, 
the term ~ ^ leads to the IR finite pieces and the term ~ Rs to the finite bremsstrahlung 
contribution. We refrain from writing down the explicit IR finite expressions for the various 
scalar form factors since they can be deduced easily from eq. (|6.6|) and the formulae spelled 
out in sections ^ and ^. 

7 Results and discussion 

7.1 Vector form factor and pion charge radius 

In fig. we show the vector form factor in the time-like region. We remark that the 
electromagnetic corrections from the reducible photon loops and from the bremsstrahlung 
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Figure 6: Vector form factor in the time-like region. The various virtual photon contributions 
are shown in comparison with the strong fourth order contribution. The latter is divided by 
a factor of 20. 



are very small, the latter contribution being so tiny that it will be discarded in what follows. 
A bit above threshold, the finite IR terms are of similar size as the irreducible photon loop 
contribution. For comparison, we have also shown the strong fourth order contribution 
divided by a factor of 20. We conclude that the electromagnetic corrections are of the 
expected size of about one percent, except for the divergence of the photon loop graphs 
at s = 4M^±. In the space-like region, shown in fig. ^, the electromagnetic effects from 
the finite IR graphs are somewhat more pronounced, the photon loop contribution is small. 
From this, we can deduce the electromagnetic corrections to the pion charge radius. 



\ / TT, strong 



\ I TT, photon 



finite IR 



31 



47r/ 5M2+ ' 



e \2 8 
4^/ iW^ 



log 



Mil 



(7.1) 
(7.2) 
(7.3) 



The numerical results for these various contributions are 



0.442 fm% (r^)^ 



strong 



photon 



-0.007 fm^ 



(r')r,finite IR = 0.024 (0.018) fm^ , (7.4) 



for = 10 (20) MeV. We use F = 93MeV, M^± = 139.57 MeV and k = 16.5. There are 
two ways of dealing with the resolution dependent IR contribution. First, if one compares 
the similar process e"'"e~ — > /i"*"/^", soft photon radiation to leading order in gives a 
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Figure 7: Vector form factor in the space-like region. The various virtual photon contri- 
butions are shown in comparison with the strong fourth order contribution. The latter is 
divided by a factor of 20. 



finite IR contribution to the "intrinsic" muon radius, which is believed to be zero 
contribution has the same form as in eq 



This 



( |7.3| ) but with the pion mass replaced by the 
muon mass (see e.g. the explicit calculation in ref. |^). If one therefore directly compares 
the processes e+e~ 7r+7r~(7) and e~^e~ — > /i+/x~(7), one is only sensitive to the difference 
(1/M^±) log(A£'/M.^±) — (1/M^) log{AE/M^) which reduces the dependence on the detector 
resolution by almost one order of magnitude. Another way of proceeding is to follow the 
arguments which have been used in the discussion of vrvr scattering | 



101. While the 



cross section of course depends on the detector resolution and so does the scattering length 
a, da/dQ = |ap, one can define an electromagnetically corrected scattering length which 
is independent of AE. This is achieved by absorbing the resolution dependent terms into 
the appropriately redefined phase space. Such a procedure applied here would essentially 
amount to dropping the term given in eq. ( [7.3|) . To summarise, the resolution dependent 
terms are probably not to be understood as "real" isospin violation effects in the reaction 
e~^e~ — > TT~^7i~, but rather as part of the radiative process e"'"e~ — *• 7r"'"7r~7, and therefore 
should not be regarded as intrinsic electromagnetic corrections to the form factor. 

We should explicitly compare these numbers to the two-loop corrections for the vector 
form factor calculated in ref. [0. Depending on the set of parameters they use, they find 
corrections at 0{q^) which numerically amount to 



2— loop 



0.017. . .0.022 fm^ 



(7.5) 



This is larger than the "intrinsic" virtual photon contribution found above by about a factor 
of 3, and thus is similar to the findings in the pion-pion scattering case, see refs. |10[] . 
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In refs. ^, also the quadratic term in the Taylor expansion of F^{s) according to 

F^{s) = l + l{r'):s + c:s' + 0{s') (7.6) 

was considered. Counterterms to fit experimental values of exactly only appear at 0{q^), 
and indeed the finding in Q is that the 0{q'^) calculation is completely insufficient to 
describe , with the final value dominated by an 0{q^) counterterm. Here, we might either 
check whether this large 0{q^) contribution can be reduced when taking into account the 
electromagnetic corrections, or we can give credit to the two-loop analysis in case these 
again turn out to be very small. The contributions found from expanding eq. ( ^.7) ) are 

= (i^^ = °-«26GeV-. (7,7) 

cLhoton = - — ) \ = -0.412 GeV-^ (7.8) 

7r,photon ^^^J 420M^± ' ^ ' 

<finiteiR = -(^) ^ log ^ = 0.807(0.595) GeV-^ (7.9) 

for AE = 10 (20) MeV, whereas, for comparison, the authors of give a value of 

= (3.85 ± 0.60) GeV"^ , (7.10) 

found by a fit of the C(g^) CHPT amplitude to the data. We see therefore that the electro- 
magnetic effects on the curvature are of a size comparable to the hadronic one-loop parts, 
although resolution dependent and independent contributions tend to cancel to some extent 
for typical values of AE. In the analysis of ref. 0, this would lead to a moderate change in 
their LEG rl^^. 

We conclude that, throughout, the intrinsic virtual photon effects on the pion vector form 
factor in the low energy region are of the expected size and smaller than the two loop effects 
calculated in refs. 0, |^. (Note finally that one particular SU{3) effect (the charged to neutral 
kaon mass difference in kaon loops) has also been found to be tiny, see ref. ||33||.) 



7.2 Scalar form factor and scalar pion radius 

We now turn to the scalar form factor. In contrast to the vector form factor, we have an 
additional electromagnetic effect, namely the pion mass difference in the hadronic loops (of 
course, there is also a tiny strong contribution to the pion mass splitting). 

Gonsider the charged pions first. The various virtual photon corrections are shown in fig. |^ 
in comparison with the strong fourth order contribution scaled down by a factor of 20. 
Here, the different contributions are somewhat more pronounced as compared to the vector 
case, but again we observe cancellations between the various effects. Note that the finite 
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Figure 8: Scalar form factor of the charged pions in the time-hke region. The various virtual 
photon contributions are shown in comparison with the strong fourth order contribution. 
The latter is divided by a factor of 20. 

bremsstrahlung contribution is again very small. The corresponding corrections to the scalar 
radius of the pion can be worked out easily (we again neglect the tiny bremsstrahlung piece) , 

- -7X^ + ' (7-12) 



('^^)^±, photon - { ..2 ' (''•l^) 



2>5 / e \ 2 8 . 



Mk<.n,..« = -(j^j ^log-g-. (7.14) 



The novel contribution due to the pion mass difference in the pion loops is denoted by the 
subscript "mass". The resolution dependent term is identical to the one obtained for the 
vector form factor. It can be dealt with in the same manner as described in the preceding 
paragraph. The corresponding values are 

.strong = 0.550 fm^ (r2)f± ^3,, = -0.0016 fm^ 
(^')^±,photon = 0.0012 fm^, (r2)f±,fi,i,, = 0.024 (0.018) fm% (7.15) 

for /S.E — 10 (20) MeV and using 1^ — 4.3. We note that the small corrections due to the 
pion mass difference (—0.3%) and the photon loops (+0.2%) nearly cancel each other to 



21 



give a tiny effect on the radius. Even the separate contributions are, again, by one order of 
magnitude smaller than the two-loop corrections ||^, 

(^')f,2-ioop = 0.016... 0.025 fm^ . (7.16) 

The normalization of the scalar form factor is also of interest. We normalize it here to the 
mass of the charged pions (although the canonical reference mass is the one of the neutral 
pion), 

f^±(0) = Ml± (1 - 0.065 - 0.012 - 0.004) , (7.17) 

where the first correction is due to the electromagnetic operator ~ C, i.e. the pion mass 
difference, the second term is the well-known hadronic shift [|I| (using /a = 2.9), and the third 
term comprises the genuine electromagnetic corrections. To arrive at this last number, we 
have used the dimensional analysis of ref. |jlO|, \k\{Mp)\ < l/(47r)^, which leads to \k'^\ < 6.3. 



Consequently, the shift due to the electromagnetic corrections is about one third of the strong 
correction to the normalization of the charged pions scalar form factor. 

As in the case of the vector form factor, we wish to give credit to the analysis of the quadratic 
term of the Taylor expansion of the form factors given in refs. The remarks made 

about the purpose of such an analysis for the one-loop electromagnetic corrections in the 
case of {s) apply equally here. From eq. ( p. 171) we find 
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e \2 Z 



- 



-0.15 GeV"^ , (7.19) 



4.,p.o.„„ = -(^)'li; = -°-13GeV-. (7.20) 
4.,,„„K = -(j^) = 0,81 (0.60) GeV- . (7.21) 

Note that the difference in the numerical value for the hadronic one-loop contribution as 
compared to ref. P occurs because they employ the charged pion mass whereas, when taking 
isospin violation into account, the correct reference mass is the neutral pion mass. The value 
extracted from experimental data is @, |^ 

4 = (10.6 ± 0.6) GeV~^ , (7.22) 

with the error taken from the range of different analyses. Again, there have to be substantial 
0{q^) contributions, the electromagnetic corrections are small and will lead only to small 
modifications in the analysis of the 0{q^) counterterms. 

We now turn to the neutral pions. Here, the only effect is due to the pion mass difference 
in the pion loops, see fig. The corresponding contribution to the radius takes the form 



/ 9^ s 12 , M^± / e \ 2 2Z , , 
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Figure 9: Scalar form factor of the neutral pions in the time-like region. The only electro- 
magnetic contribution is due to the pion mass difference. It is shown in comparison with the 
strong fourth order contribution. The latter is divided by a factor of 20. 



Numerically we have 

(^')^o,^,3s = -0.0094 fm^ , (7.24) 

using M^o = 134.97 MeV. This amounts to a 2% reduction of the hadronic result. Such 
an enhancement of the electromagnetic effects for the neutral particles as compared to the 



charged ones has also been observed in pion-nucleon scattering |15, 16, O, rsl or neutral 



pion photoproduction off nucleons Finally, we study the normalization. We find 

f^o(O) = M^o (1 - 0.013 - 0.004 - 0.001) , (7.25) 

where the first correction is due to the hadronic shift [|l|, the second one is due to the 
strong isospin breaking ~ Ij (we use 1^ = 1.0), and the third term stems from the pion 
mass difference in the loops. We note that the strong isospin breaking amounts to a 30% 
correction of the strong isospin conserving shift. 

To conclude this section, let us quote the result for the isospin-breaking one-loop contribu- 
tion to the parameter c^o, 



which is again slightly larger than for the charged pions, but small compared to the two-loop 
contribution. 
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Figure 10: Normalized form factor (s) of the energy-momentum tensor for charged 
pions. The various electromagnetic corrections are shown in comparison with the strong 
fourth order contribution. The latter is divided by a factor of 20. 



7.3 Form factors of the energy— momentum tensor 

We have calculated the various form factors of the energy-momentum tensor for charged and 
neutral pions discussed in sect. ^. We refrain from discussing all of them in detail but rather 
show two typical cases. G^^ (s) is given in fig. The electromagnetic corrections look sim- 
ilar to the ones for the scalar form factor shown in fig. |]with two notable exceptions. First, 
as noted before, the photon loop corrections diverge at s = 0, and second, the corrections 
due to the pion mass difference in the pion loops grow with increasing s after the two-pion 



threshold. The tensor form factor for charged pions is shown in fig. |11| In this case, the 
hadronic fourth order correction is simply a polynomial linear in s since there are no loop 
contributions. The size of the electromagnetic corrections is typically a few percent of the 
strong fourth order contribution, with the exception of the kinematical points where diver- 
gences appear. Finally, we note that the form factors of the trace of the energy-momentum 
tensor are very similar to the scalar form factors discussed in the previous paragraph. 



8 Summary 

In this paper, we have considered various pion form factors in the presence of virtual photons 
at next-to-leading order in SU{2) chiral perturbation theory. The pertinent results of this 
investigation can be summarized as follows: 

(1) The vector form factor of the pion has been calculated including photon loops and 
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Figure 11: Tensor form factor of the energy-momentum tensor for charged pions. The 
various electromagnetic corrections are shown in comparison with the strong fourth order 
contribution. The latter is divided by a factor of 20. 



electromagnetic counterterms. The electromagnetic corrections are typically a few 
percent of the strong fourth order contribution, cf. figs. 0. The intrinsic virtual 
photon corrections to the pion charge radius is of the order of 1%. In addition, there 
is a detector resolution dependent contribution to the radius which is of the order of 
5% for AE = 10 ... 20 MeV. We have discussed how to deal with this contribution. 

(2) We have also investigated the scalar form factor of the charged and the neutral pions. 
For the charged pions, the results are similar to the vector form factor, cf. fig. For 
neutral pions, the only effect is due to the pion mass difference in the pion loops. This 
leads to a 2% reduction of the scalar radius of the neutral pion. We have also found 
that the scalar form factors at zero momentum transfer are more strongly affected by 
isospin breaking. 

(3) The pion matrix elements of the energy-momentum tensor are parametrized by a scalar 
and a tensor form factor. We have extended the chiral Lagrangian coupled to gravity in 
the presence of virtual photons. There is only one additional electromagnetic operator 
at fourth order. Since the energy-momentum tensor can couple directly to virtual 
photons, one of the corresponding form factors shows a divergence at zero momentum 
transfer not present in the scalar form factor of the pion. Typically, electromagnetic 



corrections are small, see figs. [TO, ITT 



We conclude by noting that presently the largest uncertainty in the low energy representation 
of the pion form factors is due to the uncertainties in the low energy constants U and Iq. 
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Strong two-loop corrections or the virtual photon effects calculated here are sizeably smaller. 
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